De novo assembly and characterization of the transcriptome of Morchella esculenta growth with selenium supplementation

Although Morchella esculenta (L.) Pers. is an edible and nutritious mushroom with significant selenium (Se)-enriched potential, its biological response to selenium stimuli remains unclear. This study explored the effect of selenium on mushroom growth and the global gene expression profiles of M. esculenta. While 5 µg mL−1selenite treatment slightly promoted mycelia growth and mushroom yield, 10 µg mL−1significantly inhibited growth. Based on comparative transcriptome analysis, samples treated with 5 µg mL−1 and 10 µg mL−1 of Se contained 16,061 (452 upregulated and 15,609 downregulated) and 14,155 differentially expressed genes (DEGs; 800 upregulated and 13,355 downregulated), respectively. Moreover, DEGs were mainly enriched in the cell cycle, meiosis, aminoacyl-tRNA biosynthesis, spliceosome, protein processing in endoplasmic reticulum pathway, and mRNA surveillance pathway in both selenium-treated groups. Among these, MFS substrate transporter and aspartate aminotransferase genes potentially involved in Se metabolism and those linked to redox homeostasis were significantly upregulated, while genes involved in isoflavone biosynthesis and flavonoid metabolism were significantly downregulated. Gene expression levels increased alongside selenite treatment concentration, suggesting that high Se concentrations promoted M. esculenta detoxification. These results can be used to thoroughly explain the potential detoxification and Se enrichment processes in M. esculenta and edible fungi.


INTRODUCTION
As an essential micronutrient for humans and animals, selenium (Se) is crucial for cellular physiological processes, including cellular antioxidant defenses, regulating thyroid hormone metabolism, and cell growth.Daily selenium supplementation is therefore important for human health (Fan et al., 2024;Kieliszek, 2019;Razaghi et al., 2021).However, selenium has harmful and beneficial attributes.Selenium intake exceeding the 400 µg daily limit may cause vomiting, diarrhea, abdominal pain, fatigue, hair loss, joint pain, or increase incidence of prostate cancer or risk of type 2 diabetes (Hadrup & Ravn-Haren, 2020;MacFarquhar et al., 2010;Spiller & Pfiefer, 2007;Thompson et al., 2016).
Therefore, people must pay more attention to dietary selenium consumption to avoid selenium-deficient defects and potential overdose risks (Rayman et al., 2018;Thompson et al., 2016).
Despite their use as a dietary supplement in China and other countries, the metabolism and responses of selenium-fortified mushrooms are poorly understood.Selenite's effect on mushroom fruiting and mycelial development and its adverse impact from toxic quantities have been documented in various mushrooms (Song et al., 2022;Wang et al., 2016).For example, addition of 1.0 µg g −1 selenite in the medium resulted in up to 25% increase in mushroom yield of F. velutipes, while 50 µg g −1 selenite treatment significantly inhibit the growth of fruiting body (Dong, Xiao & Wu, 2021).Furthermore, some studies noted discrepancies in the ability of disparate mushroom species to accumulate and utilize selenite (de Oliveira, Naozuka & Landero-Figueroa, 2022;Dong, Xiao & Wu, 2021).The variations in absorption and utilization efficiency of different mushroom fruiting bodies may be attributed to the mushroom's cytoreductive system and the metabolic pathways associated with selenium metabolism (Wang et al., 2016;Zhang et al., 2022b;Zhao et al., 2020).For example, selenium uptake by F. velutipes decreased by 25 or 26% upon phosphate or sulfite starvation, respectively (Wang et al., 2016), while Se-enriched Pleurotus citrinopileatus displayed significantly higher sulfur metabolism and ABS transporter activity (Zhao et al., 2020).At the same time, there have been a number of studies on the biological processes and reduction mechanisms about Se transformation in edible mushrooms.Transcriptional profiling of Auricularia cornea under selenium treatment showed that selenium accumulation increased genes linked to glutathione metabolisms, translation, and other metabolic processes (Li et al., 2019).Furthermore, tandem mass tag-based quantitative proteomic analysis revealed that thioredoxin 1, thioredoxin reductase (NADPH), glutathione reductase, and cystathionine gamma-lyase were key proteins involved in selenite reduction and methylation in G. lucidum (Xu et al., 2023).Thus, distinct transporters and metabolic pathways may mediate and produce different forms of selenium compounds in mushrooms.However, research on these specific mechanisms requires further study.
Selenium uptake, tolerance, and metabolism pathways in M. esculenta were investigated to better understand the effect of toxic and low Se concentrations on M. esculenta and the mechanisms underlying Se accumulation.This study investigated selenium's impact on mushroom growth in M. esculenta through treatments with different amounts of selenium.Moreover, the transcriptome of M. esculenta at two specific selenite concentrations (low and high) was explored to uncover the dose-specific mechanisms of selenium enrichment and to provide a theoretical framework for future research on the metabolism and selenium enrichment of mushrooms.

Strain and culture conditions
The Morchella esculenta (L.) Pers.strain (ACCC50764) was provided by the Agricultural Culture Collection of China (Beijing, China) and was cultured on potato dextrose agar (PDA) plates (Qingdao Hopebio Technology, Qingdao, China).Solid and liquid cultures were performed in an ordinary PDA liquid medium (without sodium selenite) or Se-PDA liquid medium (supplementation with sodium selenite) at 26 ± 0.5 • C based upon experimental need.
Se tolerance and growth rate/biomass measurement Solid-state fermentation using Se-PDA or ordinary PDA plates was employed to determine growth rate and biomass production.Na 2 SeO 3 was diluted with distilled water and added to modified PDA to establish the final selenium concentrations ranging from 1 to 30 µg mL −1 .The colony diameter was measured twice, each perpendicular to the other, on the third and seventh days to evaluate the growth rate (mm d −1 ).To determine the effect of selenium treatment on the colony color, the colony color was observed against a white background.Liquid cultures were initiated by transferring three 6-mm-diam agar plugs into 15 mL of PDA liquid medium in 125 mL Erlenmeyer flasks at 26 ± 0.5 • C for two days.Next, biomass production was assessed by injecting 15 mL of the Morchella liquid strand to infect 150 mL of PDA liquid medium containing different SeO 3 2− concentrations (0.05 to 30 µg mL −1 ).The liquid culture was performed according to our previously published work (Qian et al., 2023).After 10 min of centrifugation (10,000 rpm), the mycelium was harvested and placed in an oven at 60 • C until its weight remained constant.The total Se levels within the mycelium for each treatment sample were ascertained through hydride generation-atomic fluorescence spectrometry (HG-AFS), as described by Wu et al. (2020).

Preparation of mycelial samples for transcriptome sequencing
First, 5 µg mL −1 and 10 µg mL −1 were selected as the promoting and inhibiting Se concentrations based on the selenite growth response to investigate the Se-enrichment process of M. esculenta.Following cultivation on PDA plates at 26 ± 0.5 • C for three days, the mycelial of the M. esculenta strain and the PDA media were cut into 1 cm × 1 cm pieces and placed into 250 mL flasks filled with PDA liquid medium containing 0, 5, and 10 µg mL −1 sodium selenite (labeled as CK, MSe5, and MSe10, respectively).Each sample had three replicates.After incubation at 26 ± 0.5 • C for seven days, the mycelial sample (1 g wet weight) for each treatment was harvested through filtration and stored at −80 • C until RNA extraction.

Library construction and RNA-Seq
RNA was extracted using a Total RNA Extraction Kit per the manufacturer's instructions (Solarbio Science and Technology Co., Ltd., Beijing, China).The quantify the total RNA was measured by 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA).Nine cDNA libraries were constructed using the TruSeq RNA Sample Preparation Kit (Illumina, San Diego, CA, U.S.A.) after RNA quality assessment and contaminated RNA removal.
The quality of the cDNA libraries was also measured using the 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA).The libraries were divided into three groups (CK, MSe5, and MSe10, respectively) based on the Se concentration of the treatment groups, each containing three replicates.Each cDNA library was pair-end sequenced on an Illumina NovaSeq 6000 sequencer at Shanghai BIOZERON Co., Ltd (Shanghai, China).Clean readings from raw data were obtained by filtering out ambiguous reads with more than 5% unknown bases (''N'') or duplicate reads from PCR after sequencing.

De novo transcriptome assembly and gene functional annotation
Trinity, a short-read assembly program designed for transcriptome analysis, was used for de novo assembly based on clean reads (Grabherr et al., 2011).Low-quality (Q <20) data and adapter sequences were removed from the analysis.The statistical power of this experimental design, calculated in RNASeqPower is 0.88.Next, Bowtie2 (ver.2.2.5) was used to extract mitochondrial and non-coding RNAs from trimmed pair-end reads (Langmead & Salzberg, 2012).Every cluster's longest sequence was chosen as a single unigene and sequentially annotated by comparing its transcripts with those found in Swissprot, Kyoto Encyclopedia of Genes and Genomes (KEGG) (http://www.genome.jp/kegg),Gene Ontology (GO) (http://wego.genomics.org.cn), and The National Center for Biotechnology Information's (NCBI) NR/Nt databases.Finally, GO and KEGG pathways were used to evaluate DEGs, as described by Lu & Bau (2017).

Identification of differentially expressed genes
The fragments per kilobase of exon region in a given gene per million mapped fragments approach (FRKM) was applied to calculate differentially expressed genes (DEGs) of different libraries using the DESeq R package (DESeq version: 1.30.0,R version: 4.0.4).DEGs with a significant difference were filtered using the following threshold: |log2(FoldChange)|>1 and p-value <0.05.In addition, two-way hierarchical clustering (complete linkage) was performed among the selenium-treated (MSe5 and MSe10) and control groups by adopting the R package Pheatmap to demonstrate the expression pattern of each DEG across all samples (Metsalu & Vilo, 2015).Moreover, DEGs were subjected to KEGG pathway enrichment analysis and GO functional enrichment using a previously published formula (Kim et al., 2020).A p-value cutoff of less than 0.05 was used to evaluate the significance of each TOG or KEGG pathway's enrichment pattern between the selenium-treated and control groups.

Differentially expressed gene validation using real-time quantitative PCR
Real-time quantitative PCR(RT-qPCR) was used to validate the seven randomly selected DEGs and transcriptome sequencing data reliability.Primer Premier 5.0 was used to design gene-specific primers for the chosen unigene sequences; the primer sequences are listed in Table S2.The M. esculenta strain was cultured under conditions similar to those for RNAseq library construction, as previously described.PCR reactions were performed using the one-step RT-qPCR kit (Sangon Biotech Co., Ltd., Shanghai, China) per the manufacturer's instructions.The 2 − Ct method was utilized to determine the target gene's relative expression level.The 18S rRNA was used as housekeeping gene and the expression levels of specific DEGs were normalized through comparison with the internal reference gene.

M. esculenta growth and morphology under selenium accumulation
Mycelia morphologies in selenite-supplicated media at concentrations ranging from 1 to 5 µg mL −1 were similar to those in the control group.However, visible indicators were noted with Na 2 SeO 3 concentrations over 10 µg mL −1 .The mycelia transitioned from a white to a reddish color (Fig. 1); red indicated elemental Se (Se 0 ) formation rather than organic selenium compounds in the mycelia (Xu et al., 2021a).Moreover, the growth rate of mycelium and M. esculenta biomass under different selenium treatments varied significantly by Se concentration (Table 1).The initially augmented mycelia growth was inhibited by an increase in Na 2 SeO 3 concentration.For example, adding small amounts of selenite (5 µg mL −1 ) to the medium slightly promoted mycelia growth, resulting in a 4.93% increase in mushroom yield (Table 1).However, 10 µg mL −1 selenite treatments significantly reduced mycelial development and biomass compared with the control group.These results corroborate those reported in other studies (Kim et al., 2014;Song et al., 2009;Xu et al., 2021a).For example, lower Se concentrations of 1, 10 and 100 µM significantly improved the biomass levels of the fungus Pleurotus eryngii, while higher Se concentrations of 1,000 and 10,000 µM can significantly inhibit the mycelial growth (Kim et al., 2014).Furthermore, Se concentrations of 0-80 µg mL −1 promoted mycelial growth of Stropharia rugoso-annulata, whereas the mycelia biomass decreased rapidly when the concentration exceeded 200 µg mL −1 (Song et al., 2009).Different selenium concentrations may influence mushroom growth due to the toxic effects caused by selenium exposure.Se exposure can induce toxicity and alter the growth of organisms by causing lipid peroxidation, producing reactive oxygen species (ROS), and diminishing sulfur metabolism (Kieliszek et al., 2020;Li et al., 2007;Staneviciene et al., 2023).Thus, the mushroom's tolerance and capacity to metabolize selenium may be related to its species and antioxidant system.At the same time, the life cycles and substrate composition of different edible mushrooms may also account for the differences in selenite tolerance.

M. esculenta transcriptome sequencing and de novo assembly
Gene expression profiles of M. esculenta were compared between 5 µg mL −1 (MSe5) and 10 µg mL −1 (MSe10) selenium treatments and the control (CK) using high-throughput RNA sequencing based on the Illumina HiSeq4000 platform.Since M. esculenta's genomic data has yet to be reported, de novo assembly was selected for this study.An average of 2.76×10 7 clean readings per sample were obtained after quality examination and filtration (Table S1).Thereafter, 73,135 transcripts and 73,135 unigenes were created by assembling the clean readings.Unigene lengths ranged from 227 to over 46,000 bp, averaging 2,819 bp.The number and quality score of aligned reads were generally sufficient for subsequent transcriptome analysis.

Transcriptional data analysis
First, 73,135 transcripts from the Trinity assembly were screened for subsequent DEG identification.For samples treated with 5 µg mL −1 of selenium (MSe5), 452 upregulated genes and 15,609 downregulated genes were identified (Fig. 2A).Following 10 µg mL −1 selenium treatment (MSe10), 14,155 genes were statistically upregulated (800 DEGs) or downregulated (13,355 DEGs) relative to the controls (Fig. 2B).Furthermore, co-expressed DEGs were identified using a Venn diagram and the results revealed that 475 and 2,381 genes were specifically expressed in the two comparison groups (MSe10 vs. CK, MSe5 vs. CK) and 13,680 DEGs were shared (Fig. 3).The relative expression of seven randomly selected DEGs was determined through RT-qPCR to verify the RNA-seq results (Table S2).Overall, regulation of the tested genes was consistent with the transcriptome analysis results, indicating that the latter was reliable.GO enrichment analyses demonstrated that these DEGs could be classified into three categories: cellular components (CC), biological processes (BP), and molecular functions (MF).In the MSe5 vs. CK group, the upregulated genes of cellular components category were primarily part of the membrane's cellular components (111 DEGs) and cell parts (93 DEGs).Regarding biological processes, upregulated genes were implicated in metabolic processes (85 DEGs), cellular processes (85 DEGs), and biological regulation (25 DEGs).Furthermore, upregulated genes classified into the molecular function category were associated with binding (143 DEGs), catalytic activity (120 DEGs), and transporter activity (24 DEGs) (Fig. 4A, Table S1).
The GO functional annotation results indicated that selenium enrichment or treatment can regulate numerous genes related to energy production/metabolism, protein synthesis, cell wall organization or biogenesis, and carbohydrate catabolism (Jiao et al., 2022;Zhou et al., 2018).Notable, the mostly highly enriched GO terms between the MSe5 vs. CK and MSe10 vs. CK groups were identical, suggesting that genes in these clusters were clearly Se-related.For example, DEGs participated in biological regulation or binding were all highly enriched in the MSe10 vs CK group and MSe5 vs. CK group.Similar to our results, Li et al. (2023) found that feeding selenium-enriched yeast-treated diets can significantly regulate the expression of genes involved in metabolism-related biological processes, such as biological regulation, molecular processes in Broiler chickens (Gallus gallus).While Se treatment can modulate the global transcriptome of Yorkshire-Landrace gilts.The significantly enriched GO terms in selenite feeding group are mainly related with protein binding and translation initiation factor activity (especially translation factor activity-RNA binding) (Dalto et al., 2018).Meanwhile, higher selenium treatments upregulated more genes than lower concentration treatments.Other studies reported similar findings regarding selenium enrichment in Cardamine hupingshanensis (Zhou et al., 2018), and Bombyx mori (Jiang et al., 2020).For example, in the hyperaccumulator plant C. hupingshanensis, more DEGs (42 upregulated and 61 downregulated) with various physiological processes and cellular functions were apparent following high Se concentration treatment (80,000 µg Se L −1 ).Contrastingly, only 23 upregulated and 7 downregulated genes were found in the low Se (100 µg Se L −1 ) treatment group (Zhou et al., 2018).
Quantitative proteomic analysis of Ganoderma lucidum cultured with 200 ppm selenite showed that a number of metabolic pathways were prominently enriched such as pyruvate metabolism, citrate cycle, tyrosine metabolism, and glycolysis/gluconeogenesis (Xu et al., 2023).Therefore, KEGG enrichment analysis indicated that these genes may play roles on Se accumulation or tolerance in M. esculenta.It is noteworthy that several DEGenriched pathways in the MSe5 vs. CK and MSe10 vs. CK groups were identical, primarily encompassing cell growth and death (21 DEGs), translation (17 DEGs), aging (five DEGs), and carbohydrate metabolism (four DEGs).Similar results were also observed following selenium accumulation in the Camellia sinensis tea plant.Under selenite treatment, DEGs mainly enriched in ribosomes (52 unigenes) and oxidative phosphorylation (14 unigenes) were upregulated in the roots and leaves (Cao et al., 2018).In the fat body of silkworms after 50 µM and 200 µM selenite treatment, DEGs were predominantly associated with the KEGG pathways, including hypertrophic cardiomyopathy, dilated cardiomyopathy, and the VEGF signaling pathway, and salivary secretion (Jiang et al., 2020).The consistent expression of genes under different treatment conditions indicates that selenite treatment activated genes in these pathways.

Genes involved in Se enrichment or biotransformation
Uptake is the first step for cells to assimilate and utilize selenium; thus, efficient transport systems are crucial for selenium uptake (Wang, Rensing & Zheng, 2022a).Several transporters involved in Se uptake and translocation have been identified (Wang et al., 2024).The monocarboxylate transporter Jen1p and phosphate transporters Pho87P/Pho90P/Pho91P, which facilitated selenite with high affinities, are responsible for selenite accumulation in yeast (McDermott, Rosen & Liu, 2010).In addition, selenite can be transported by several bacteria, including Escherichia coli and Proteus sp.YS02., via the sulfate ABC transporter complex (Rosen & Liu, 2009;Wang et al., 2022b).On the other hand, selenite transporters within edible fungi have not been studied at the molecular level.This study found no significant changes in gene expression, such as sulfate transporters or phosphate transporters, but revealed that the MFS substrate transporter (NODE_1033) was significantly upregulated following selenite treatment (MSe5 vs. CK: 5.15 Log 2 FC, MSe10 vs. CK: 5.66 Log 2 FC).The upregulation of MFS transporter gene expressions has been reported in selenium-enriched lactic acid bacteria (Liao & Wang, 2022).Thus, the continuous enhancement of MFS substrate transporter (NODE_1033) gene expressions in this study indicated that non-specific transporter proteins may carry Se oxyanions in the culture medium into the M. esculenta cells.Moreover, gene expressions related to redox homeostasis, such as genes encoding thioredoxin and glutathione S-transferase (GstA), peroxisome, and glutathione synthetase, were significantly expressed following treatment with different selenite concentrations (Jiang et al., 2020;Wang et al., 2023).For example, the Log 2 FC of genes encoding thioredoxin and glutathione S-transferase were 10.2 and 9.27 in the MSe5 vs. CK group, respectively.In the MSe10 vs. CK group, the Log 2 FC of these gene expressions increased to 10.92 and 10.49, respectively.These findings indicate that as selenite treatment concentrations increase, the gene expressions are also upregulated.In comparison, genes encoding glutamate-cysteine ligase (GCS) were only significantly expressed in the MSe10 vs. CK group (Log 2 FC = 7.16).
Reactive oxygen species (ROS) can damage cell membranes (Lu et al., 2022).Thioredoxin (TRX) and glutathione are integral for ROS detoxification during stress responses and for determining the thiol/disulfide status of proteins by sensing and reducing equivalents to many target proteins, including reductases and peroxidases (Zhou et al., 2018).It is commonly known that oxidative stress predicates selenium cytotoxicity.If the selenium concentration in the environment is too high, selenium stress can activate the antioxidant system to scavenge ROS, enhancing the resistance of microbes or plants to oxidative stress (Luo et al., 2021;Varlamova et al., 2021).Therefore, increasing the expression of these genes demonstrated that these proteins are crucial for maintaining the redox balance in M. esculenta cells under selenium stress, establishing cellular resistance to selenium.
Mushrooms absorb selenium from the soil or culture medium and convert it into organic Se (primarily selenoproteins and selenopolysaccharides) (de Oliveira, Naozuka & Landero-Figueroa, 2022;Qian et al., 2023).Se methylation and selenocysteine (SeCys) generation are the two main stages of selenoprotein synthesis.Selenite treatment upregulated gene expressions involved in cysteine and methionine metabolisms, whereas those enriched in cysteine synthesis pathways were significantly downregulated.For example, the expression of genes encoding aspartate aminotransferase, which possess cysteine sulfinate and aspartate aminotransferase activities, were significantly upregulated following 10 µg mL −1 selenite treatment with a Log 2 FC value of 5.5.This result is consistent with previous findings concerning selenium accumulation in other eukaryotic organisms (Cao et al., 2018;Zhang et al., 2019).For instance, genes involved in the metabolism of sulfur, cysteine, and methionine were markedly upregulated in the roots of the tea plant Camellia sinensis after selenite treatment and are probably linked to Se accumulation (Cao et al., 2018).Under acid stress, genes related to ATP production and sulfur metabolism were significantly upregulated, including methionine, cysteine, and GSH production in selenium-enriched Candida utilis (Zhang et al., 2017).
Previous studies have discovered that the total flavonoid and chlorogenic acid contents of several high plants can be increased if the applied dose of selenium does not exceed a certain concentration (Guo et al., 2020;Li et al., 2020).Therefore, the expression profiles of structural genes associated with isoflavones and polysaccharides were examined to determine whether sodium selenite treatment would enhance the biosynthesis of either compound in M. esculenta.The results indicated that two genes related to isoflavone biosynthesis or regulation, eight related to flavonoid metabolism, and five in biosynthetic pathways were downregulated (See Table S5 and Table S6).For example, the expression of genes encoding phosphoglucose isomerase, which directs Glc-6-phosphate conversion, was significantly downregulated following 5 µg mL −1 and 10 µg mL −1 selenite treatment (Log 2 FC = -13.53for both).Furthermore, the expression of genes encoding flavanone 4-reductase/Bifunctional dihydroflavonol 4-reductase, involved in flavonoid metabolism, was significantly downregulated at 10 µg mL −1 selenite treatment (Log 2 FC = −9.6)but showed no effect following 5 µg mL −1 selenite treatment.These results are consistent with previous study indicating that an application dose exceeding 2.0 mg kg −1 impairs total flavonoid and chlorogenic acid levels (Wang et al., 2021).Furthermore, an initial rise of total polyphenol contents/total flavonoid content followed by a reducing trend was detected in selenium biofortified germinated black soybeans with the increases in Se concentration (Huang et al., 2022).However, contrary to our results, external sodium selenite can significantly induce the expressions of genes involved in flavone and flavonol biosynthesis (Li et al., 2020).The current study's results deviate from previous research and emphasize the significance of employing proper selenite concentrations for treatment.Appropriate sodium selenite applications can increase the levels of selenium and other active chemicals in eukaryotic organisms; however, it may have the opposite effect after surpassing a specific concentration threshold.

CONCLUSIONS
In the food industry, selenium-enriched mushrooms play a significant role as a daily supplement.Therefore, RNA-Seq was used to evaluated the effect of different concentrations (5 µg mL −1 and 10 µg mL −1 ) of Se on the gene expression changes of M. esculenta.Both treatments can cause the expression level changes of genes involved in cell cycle, meiosis, aminoacyl-tRNA biosynthesis, spliceosome, mRNA surveillance pathway, and protein processing in endoplasmic reticulum pathway.Among which, one MFS substrate transporter and one aspartate aminotransferase genes were significantly upregulated and were presumed to be involved in Se metabolism.Moreover, the level of gene expression changed increased with selenite treatment concentration, suggesting that the high levels of Se promoted the detoxification of M. esculenta.However, genes related to isoflavone biosynthesis and flavonoid metabolism were significantly downregulated under selenite stress, highlighting the importance for selenite treatment with appropriate concentrations to obtain the balance of Se content and other active compounds.The results may serve as a basis for elucidating the specific effects of Se on growth and accumulation/detoxication processes in M. esculenta and promoted functional seleniumenriched mushroom production.

Figure 1
Figure 1 The morphologies of mycelia of M. esculenta grown on media containing different concentrations of selenite.The morphologies of mycelia of M. esculenta grown on media containing selenite ranging from 1 to 30 µg mL −1 .

Figure 2
Figure 2 Volcano plot of DEGs between (A) MSe5 vs. CK and (B) MSe10 vs. CK groups.The diamond, circle and square pots indicate differentially upregulated, differently downregulated, and non-differentially expressed genes, respectively.Each point in the diagram represents a gene.Full-size DOI: 10.7717/peerj.17426/fig-2

Figure 4
Figure 4 GO functional classification of M. esculenta unigenes different groups.GO functional classification of M. esculenta unigenes in (A) Upregulated genes in MSe5 vs. CK, (B) Downregulated genes in MSe5 vs. CK, (C) Upregulated genes in MSe10 vs. CK, and (D) Downregulated genes in MSe10 vs. CK groups.The X -axis depicts the number of DEGs, and the Y -axis correlates to various gene functions.CC, cellular components; BP, biological processes; MF, molecular functions.Full-size DOI: 10.7717/peerj.17426/fig-4

Figure 5
Figure 5 The KEGG enrichment analysis of DEGs within different M. esculenta groups based on enrichment factors.(A) The top nine pathways of upregulated genes in MSe5 vs. CK; (B) The top 13 pathways of downregulated genes in MSe5 vs. CK; (C) The top 13 pathways of upregulated genes in MSe10 vs. CK; (D) The top 15 pathways of downregulated genes in MSe10 vs. CK.The horizontal axis represents the name of the enriched KEGG pathway, the vertical axis represents the enrichment factor of significantly different genes in that pathway, the dot color represents the Q-value, and the dot size represents the number of genes.Full-size DOI: 10.7717/peerj.17426/fig-5

Table 1 The effect of sodium selenium on the growth and morphology of M. esculenta.
*Mean ± standard error.